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ABSTRACT 

A method for the asteroseismic analysis of f3 Cephei stars is presented and 
apphed to the star u Eridani. The method is based on the analysis of rotational 
splittings, and their asymmetries using differentially-rotating asteroseismic mod- 
els. Models with masses around 7.13 Mq, and ages around 14.9 Myr, were found 
to fit better 10 of the 14 observed frequencies, which were identified as the fun- 
damental radial mode and the three i = 1 triplets gi, pi, and p2. The splittings 
and aymmetries found for these modes recover those provided in the literature, 
except for p2. For this last mode, all its non-axysimmetric components are pre- 
dicted by the models. Moreover, opposite signs of the observed and predicted 
splitting asymmetries are found. If identification is confirmed, this can be a very 
interesting source of information about the internal rotation profile, in particular 
in the outer regions of the star. 

In general, the seismic models which include a description for shellular ro- 
tation yield slightly better results as compared with those given by uniformly- 
rotating models. Furthermore, we show that asymmetries are quite dependent 
on the overshooting of the convective core, which make the present technique 
suitable for testing the theories describing the angular momentum redistribution 
and chemical mixing due to rotationally-induced turbulence. 

Subject headings: stars: evolution -stars: individual: u Eridani — stars: in- 
teriors - stars: oscillations (including pulsations) stars: rotation -(Stars: vari- 
ables:) P Cephei 
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Introduction 



The (3 Cephei star u Eridani (HD 29248) is nowadays one of ttie most in-deptli studied 
stars. Classified as a B2III star, it presents a relatively simple internal structure, char- 
acterised by a large conve ctive core. The k mechanisrn located in the metal opacit y bump 
situated around 2 x 10^ K (jPziembowski &: Pamiatnykhlll993l : iGautschy fc Saidll993l ). drives 
its oscillations. In addition, the pulsation periods of /? Cephei stars, varying in the range of 
3-8 hours, makes them suitable for detecting and analysing oscillation frequencies. Several 
photometric and spectroscopic multisite observational campaigns for i / Eridani have bee n set 
up since 2002 wit h subsequent frequency analysis from spectrosc opy (lAerts et al.ll2004l ) and 
from photometry (IHandler et al.l 12004 : iJerzykiewicz. et al.ll2005l ). In these later works, two 
independent low-frequency, high-order g modes were detected. As a consequence, u Eridani 
can also be classified as a SPB star. In terms of light curves, these multisite observations 
constitute the largest time-series ever collected for a (3 Cephei star. Nowadays, u Eridani 
presents the richest oscillation spectrum (14 independent frequencies) of the (3 Cephei-type 
class. Such privileged scenario for asteroseismology has led to several authors to analyse its 
oscillation spectrum and perform seismic models of the star. 

Several attempts have been carried out to provide a plausible sei smic model which ex- 
plains the observed frequencies of u Eridani. In lAusseloos et al.l (120041 ) (hereafter ASTA04), 
a massive exploration of standard and non-standard stellar models was undertaken in order 
to fit the oscillation data. The authors showed that an increase in the relative number frac- 
tion of iron throughout the whole star, or a large decrease in the initial hydrogen abundance, 
made the stellar models satisfy all the observational constraints, in particular, the modes 
around the fundamental radial mode are predicted unstable. 



Pamyatnykh et al.l (120041 ) performed a seismic analysis of the oscillation spectrum of 
u Eridani taking the excitation of modes into account. In that work only three frequencies 
were fitted, failing to reproduce mode excitation in the broad observed frequency range of the 
{i = 1,^2) modes, associated with the highest frequency peak in the spectrum. Nevertheless, 
they also inferred some properties of the internal rotation rate using the rotational splittings 
of two dipole (£ = 1) modes identified as gi and pi. In particular, their results suggest that 
the mean rotation rate in the core and the /i-gradient zone is about three times higher than 
in the envelope, for their two standard models fitting the three aforementioned frequencies. 



Recently, jPziembowski &: Pamyatnykh] (120081 ) have analysed the impact of considering 
uncertainties in the opacity and element distributions on the interpretation of u Eridani's 
oscillation spectrum. No satisfactory explanation of the low-frequency modes was found. 
Moreover, the authors concluded that some enhancement of the opacity in the driving zone 
is required. 
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The rotational splitti ng; asymmetries of v Eridani have al so been studied under different 
hypothesis. In particular. iDziembowski &: Jerzykiewicz J2003h suggested that the asymmetry 
of the i = \ triplet (around 5.64 d~^), as measured by lvan Hoof. A.l(ll96ll ). could be explained 
by two principal effects: the quadratic effects of rotation and a s trong magneti c dipo le field 
of the order of 5-10 kG. Such a magnetic field was searched for by lSchnerr et al.l (120061 ) using 
spectropolarimetry with no success. 
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Table 1. A list of independent frequencies taken from I Jerzykiewicz. et al.l (120051 ). 



ID 




fi [/uHz] 


/i 


5.7632828 


66.7047 


/2 


5.6538767 


65.4384 


/3 


5.6200186 


65.0465 


h 


5.6372470 


65.2459 


h 


7.898200 


91.4143 


h 


6.243847 


72.2667 


/r 


6.262917 


72.4875 


h 


7.20090 


83.3437 


h 


7.91383 


91.5952 


fio 


7.92992 


91.7815 


hi 


6.73223 


77.9193 


fl2 


6.22360 


72.0324 
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Motivated by these results, the present work aims at performing a complete modelling 
of V Eridani taking the effect of rotation up to second order into account, with the special 
feature of considering the presence of a radial differential rotation in the seismic modelling. 
To do so, a method based on the analysis of rotational splittings and their asymmetries is 
discussed. 

The paper is organised as follows: Section [2] describes the modelling procedure and 
provide details of both the evolutionary models and the oscillation spectra computation. In 
Section [31 the different sources for constraining the stellar parameters are compared, which 
include a stability analysis. Then, Sections H] and [5] explain the method here presented and 
discuss its application to the particular case of v Eridani. Finally, conclusions and final 
remarks are written in Section O 



2. Seismic modelling 

The seismic modelling described below consists in the computation of evolutionary mod- 
els and their corresponding adiabatic and non-adiabatic oscillation spectra. This is described 
in the following sections. 



2.1. Equilibrium models 



To theoretically characterise v Eridani, we build e quilibrium models representative of 
the star with the evolutionary code cesam (jMorellll997l ). In particular, models taking first- 
order effects of rotation into account are constructed. Such models are the so-called pseudo- 
rotating models, whose spherically averaged contribution of the centrifugal acceleration is 
included by means of an effective gravity gcg = g — Ac{r), where g is the local gravity, r is 
the radius, and Ac{r) = 2/3 rfi^(r) is the centrifugal acceleration of matter elements. This 
spherically averaged component of the c entrifugal acceleration does ri ot change the order 
of the hydrostatic equilibrium equations ( Kippenhahn fc WeigertlEoOol ). The non-spherical 
components of the centrifugal acceleration (which are not included in the equilibrium mod- 
els), are included in the adiabatic oscillation computations (see next section) by means of 
a linear perturbation analysis according to Soufi, Goupil & Dziembowski (1998) (see also 
Suarez, Goupil & Morel 2006). It is possible to evaluate the impact o f a differential rota tion, 
using two simple hypothesis when prescribing the rotation profile (jSuarez et al. 2006 ): 1) 
instantaneous transport of angular momentum in the whole star (global conservation) which 
thus yields a uniform rotation, or 2) local conservation of the angular momentum (shellular 
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rotation), except in the convective core whose rotation is assumed to be rigid. In both cases, 
no mass loss is considered at any evolutionary stage, that is, the total angular momentum 
is assumed to be conserved. These hypothesis represent extreme cases, so reality is presum- 
ably somewhere in between. In fact, similar rotation profiles have been found when analysing 
the evolution of giant stars inclu ding rotationally induced mixing of chemical elements and 
transport of angular momentum (IMaeder fc Meynetll2004l ). 



Input physics have been adequately chosen for main sequence B star s. Particularly, for 

themass range treated in this work (7-13 Mq), the CEFF equation of st ate (IChristensen-Dalsgaard &: Daep] 



1992 ) is used, in which, the Coulombian correction to the classical EFF ( Eggleton et al. 1973h 



has b een included. The opacity tables are taken from t he OPAL package Jlgles ias fc Rogers 



19961 ). For the metal mixture, the abundances given bv lCrevesse fc Noe ls (199$" are used 



A weak electronic s creening is ass umed, which is valid in the evolutionary stages con- 



sidered in this work (see IClaytonlll968l . for more details). For adiabatic stellar models, the 



Eddington's T{t) law (grey approximation) is considered for the atmosphere reconstruction. 
For non- adiabatic stellar mo dels, the atmo spheres are reconstructed using the Kurucz equi- 
librium atmospheric models (lKurucall993[ ) from a specific Rosseland optical depth until the 
the last edge (around r = 10~^) of the star is reached. 



Convection is treated with the mixing-length theory (iBohm-Vitensd Il958l ) which is 
parametrised with omlt = ^/Hp, where / is the mean path length of the convective elements, 
and Hp is the pressure scale height. In addition, we use the overshoot parameter, defined as 
dov (^ov being the penetration length of the convective elements). 

The grid of equilibrium models has been constructed with steps of 0.002 dex, 0.01, and 
0.05, in the metallicity Z, overshoot parameter dov, and mass M, respectively. 
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Fig. 1. — Mass-Metallicity relation fitting the observed frequencies /i and /4 obtained from 
three different codes: CLES, GRANADA PACKAGE and warsaw-new jersey codes. Two 
comparisons are shown: 1) for models without overshooting (liege package, GRANADA 
PACKAGE & WARSAW-NEW JERSEY) and 2) for models with dov =0.2 and 0.3 (liege pack- 
age, GRANADA PACKAGE). 
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2.2. Oscillations 



Adiabatic eigenfrequencies of selected pseudo-rotating models described in the previous 



section have been computed with the oscillation code FILOU (see iTran Minh fc Leonlll995 



Suareal2002l ). In this code, oscillation frequencies are obtained by means of a perturbative 
method taking into account up to second order effects of rotation. In the case of near- 
degenerate frequencies, i.e. when two or more frequencies are close to each other ( u jnim. y 
^^n'l'm), the corrections for near degeneracy are included. As detailed in lSuarez et al.l (120061 ) . 
the oscillations computation also takes into account the presence of a radial differential 
rotation profile under the form 

Q{r)=Q{l + r]o{r)^ (1) 

where Cl represents the angular rotational velocity at the surface and i]Q{r) a radial function. 
This rotation profile is equivalent to the shellular rotation profile obtained with the pseudo- 
rotating models described in the previous section. 

Concerning the insta bility computation s, non-adiabatic theoretical observables are ob- 
tained using GRACO code (IMoya et al.ll2004j ). This code solves the stellar pulsation equation 
in a non-adiabatic non-rotating frame by dividin g the star in two p arts, 1) the interior, by 
means of the non-adiabatic equations described in lUnno et al.l ( 119891 ). a nd 2) the atmosphere , 
taking into account the interaction with the pulsation as prescribed by lDupret et al.l (120021 ). 



2.3. Comparison of numerical seismic packages 



As it is widely known, in stellar modelling, and in particular, in the field of stellar seis- 
mology, the use of different codes with di fferent numerical techniques, can be crucial for the 
correct interpretation of seismic data (see lMoya et al.ll2008l ). ASTA04 compa red two seismic 



packa ges: the liege pagkage, which co mprises the evoluti onary code CLES ( IScuflaire et al. 
20071 ). plus the os cillation code LOSC (iBourv et al.l Il975r) . and the warsaw-new JER- 



SEY PAGKAGE (by iDziembowski fc Jerzykiewicj l2003l ) which comprises the warsaw-new 
JERSEY code and its corresponding oscillation code. Following ASTA04, these codes are 
compared (through a mass-metallicity relation) with our seismic model ling package , the 
GRANADA PAGKAGE, which is composed by th e evolutionar y code gesam (|Morellll997r) . and 



the oscillation code GRAGO (IMoya et al.ll2007l ) and FiLOU (ISuarez et al.l 120071 ) . described in 
Sections 12.11 and 12.21 respectively. Hereafter, those three packages are called LP, WP, and 
GP, respectively. 

Similarly as done in ASTA04, the mass M, metallicity Z, overshoot parameter (iov, and 
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initial hydrogen content Xi, are then varied such as to fit the two observed frequencies, /i 
and /4 (Table [1]) by using non-rotating models. This procedure yields a mass-met allicity 
relation for each ciov value, which is compatible with the similar study described in ASTA04 
(see Fig. [T]). It is found that, under similar conditions, i.e. for a given metaUicity, q;mlt 
and dov, GP predicts higher mass values, around 0.25 M©. These differences increase slightly 
when dov increases. Such a different behaviour can be explained by the different treatments 
of the overshooting implemented in the evolutionary codes. Indeed, for the overshooting de- 
scription, the Liege evolutionary code only takes the density variations (in the overshooted 
region) into account, which slightly affects the temperature gradient, whereas cesam consid- 
ers an additional restriction by imposing that the real temperature gradient must be equal 
to the adiabatic one, i.e. V = Vad- This would imply that, either transport of heat is purely 
radiative, or it is efficiently transported outwards from the stellar core through convective 
movements (as a result of the overshooting), respectively. This constitutes an interesting 
challenge for asteroseismology because the oscillation modes are sensitive to the physical 
description of the /x-gradient zone. 
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4.40 4.39 4.38 4.37 4.36 4.35 4.34 4.33 

log T^fj. 

Fig. 2. — HR diagram containing the two estimates for the effective temperature and lumi- 
nosity (filled star symbols) adopted in this work for v Eridani (more details in Section [3?T|) . 
with the corresponding photometric uncertainty boxes (shaded area). For illustration, three 
evolutionary tracks of 7.60, 7.13 and 7M0, are also depicted. Those models have been com- 
puted using an initial hydrogen content of Xi = 0.5. The filled circle represents the selected 
model NRl, which predicts unstable and better fits the observed frequencies /i, f^, /e, and 
/g. More details Section [3l2l 
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3. Constraining stellar parameters 

The modelling of any star entails the constraint of its stellar parameters. For doing so, 
the choice of the set of free parameters to be fixed generally depends on the observational 
material available. In the case of v Eridani, the following space of free parameters was 
chosen 

V = r{M,t,do.,Z;n) (2) 

where M is the stellar mass, t the age, dov the overshooting parameter, Z the metallicity, 
and fl the angular rotational velocity. From photometric and spectroscopic observations, we 
search for an estimate of the star location in the HR diagram. This allows us to constrain the 
metallicity and to have an estimate of the evolutionary stage of the star. Stability analysis 
significantly reduces the region of the HR diagram in which representative models can be 
searched. 

Then, the fitting of four of the observed frequencies permits to better constrain the 
mass, metallicity, evolutionary stage (age) and overshooting parameter of the star. The 
small observed vsinz of u Eridani allows the use of non-rotating models for this exercise. 
However, further model constraining makes it necessary to take the stellar rotation into 
account. In particular, seismic models including shellular rotation profiles are used (see 
Section [2]). Finally, analysis of the rotationally-split modes and their asymmetries allows to 
make a refined search for representative models of the star. 



3.1. Locating u Eridani in the HR diagram 



In ord er to locate the s t ar in the HR diagram we followed the works by iMorel et al. 

(120061 ) and iDe Ridder et al.l (120041 ). based on high-precision spectroscopy. The error box 
shown in Fig. [2] takes into account the results reported in both papers. 
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Fig. 3. — Growth rates as a function of the oscillation frequency. Each curve represents the 
growth rate obtained for oscillation modes with different degree {i ranging from to 3) . The 
lowest frequencies correspond to the SPB-type pulsation frequencies. Only the frequencies 
around the fundamental radial mode are predicted unstable (positive growth rate), which 
are those kept for the present investigation. 
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On the other hand, the location of stars in the HR diagram depends on both the 
rotational velocity of the star and the inclination angle with respect to the observer, i. 
A technique to estimate and correct for the effect of fast rot ation on the deterrn ination 
of fundamental parameters for pulsating stars is describ ed in iMichel et al.l (119991 ) . This 
technique was then refined by lPerez Hernandez et al.l (119991 ) . who applied it to S Scuti stars in 
clusters. A first consequence that can be extracted from those works is that rotation increases 
the size of the uncertainty box of stars in the HR diagram (when obtained from photometry). 
In general, such additional uncertainties increase for increasing rotational velocities (for a 
given angle of inclination of the star). For u Eridani, a projected rotation al velocity of 
ysim_~ 16kms~^ has been observed. This value has been derived by iMorel et al 



( 120061 ) from high-resolution spectroscopy, taking into account the broadening 
due to oscillations. In principle, it can be considered as a slow rotator so that the 
additional uncertainties coming from the effect of rotation can be considered within the 
already large photometric uncertainties. 
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Table 2. Characteristics of the best models selected in this work. The different columns 

are, from left to right: the model identification, the metaUicity Z, the mass M, the radius 
R, the logarithm of the effective temperature Tcs (in K), the logarithm of the gravity g (in 

cgs), the logarithm of the luminosity (cgs), the overshooting parameter (in the Hp scale), 
the rotational frequency of the surface Us (in fiHz), the rotational frequency of the core Uc 

(in /xHz), the age (in Myr), the hydrogen abundance in the core Xc, and the radius of the 

convective core Vc- 



Model 


Z 


M/Mq 


R/Rq 




logy 


log(L/L0) 


dov 


Us 




Age 




rc/R 


NRl 


0.019 


7.13 


5.704 


4.364 


3.778 


3.923 


0.28 








14.82 


0.139 


0.129 


SRI 


0.019 


7.13 


5.714 


4.364 


3.777 


3.923 


0.28 


0.258 


0.681 


14.90 


0.138 


0.129 


SR2 


0.019 


7.13 


5.727 


4.362 


3.775 


3.918 


0.24 


0.301 


0.778 


14.60 


0.129 


0.125 


URl 


0.019 


7.13 


5.511 


4.368 


3.808 


3.909 


0.28 


0.377 


0.377 


15.80 


0.159 


0.136 
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3.2. Instability predictions 

Once the star is located in the HR diagram, models which predict the observed fre- 
quencies to be unstable (considering the constraints on physical parameters given above) are 
searched for. 

To do so, no specific ad-hoc modification of the iron mixture throughout the star is 
proposed. The search for the best models is then enhanced by using the M — Z relations 
presented in Section 12. 3[ We recall that lines in Fig. [1] represent the models which fit 
simultaneously both /i and f^. These models constitute a first guess to the best solution. 
Finally, our best solution is given by those models that also fit /e and /g. Consistently 
with ASTA04, when using standard solar mixture for the metallicity, an initial hydrogen 
content of ~ 0.50 is required in order to predict unstable the observed frequencies. In 
Fig. [21 the different evolutionary tracks shown have been computed using Xi = 0.50 in the 
ZAMS. From our grid of models (see Section [2.11) . the best model (filled circle in Fig. [2]) 
was selected to be the one which predicts unstable and better fits the observed frequencies 
/i) /4) fe, and /g. This best non-rotating model is characterised by a mass of 7.13 Mq, a 
solar metallic mixture (with an initial hydrogen content of Xi = 0.50), and by the physical 
parameters aMLT = 1, dov = 0.28 (NRl in Table [2]). Note that rotation has not been taken 
into account in the instability predictions described above. The lack of theories describing 
the effect of rotation on mo de stability rnakes it difficult to estimate this effect for u Eridani. 



Nevertheless, as stated by [Pamyatnykh[ ([19751). mode stability depends predominantly on 



the effective temperature of the models (see [Suarez et al.[[2007[ . for an interesting discussion 



on this issue), whose variations due to rotation are expected to be small (due to the low 
rotational velocity of the star). 

The results of the stability analysis obtained for that model are depicted in Fig. [3l in 
terms of growth rate as a function of the oscillation frequency. As expected, in the frequency 
region around the fundamental radial mode, the results obtained for the different spherical 
degrees i are quite similar, which implies that this parameter cannot be discriminated. 



4. Procedure 

As discussed in Sections 12. 1112. 2[ the small rotational velocity of the star makes it plau- 
sible to initially adopt the parameters of the non-rotating models which were found to be 
representative of the star. In particular, in order to work with models located in 
the HR error box which predict unstable the observed frequencies, a value of 
Xi = 0.5 for the initial hydrogen fraction is kept (see previous section). Indeed, 



- 16 - 



such a value is rather unrealistic. This can be solved ei ther by considering 



an ad- 



hoc iron enhancement in the driving zone (as done by iPamyatnykh et al.ll2004l ). 
or by modifying (uniformly throughout the stellar interior) the relative number 
fraction of iron (as done by ASTA04). However, the modelling techniques used 
here do not allow to make such modifications. Instead, following ASTA04 we 
simulate the metallicity change by modifying the initial hydrogen fraction. This 
yields, similarly to ASTA04 (see Table 1 in that paper), models around 7Mq and 
Z = 0.018 — 0.019. Similarly to previous works, standard models, i.e. those with 
Xi ~ 0.7, within the HR error box present masses higher than 8M0, Z ~ 0.015, 
for similar surface rotational velocities. Such variations with respect to the non- 
standard models have no significant impact neither on the rotation profile nor 
on the splitting asymmetries. 

A systematic search for representative models within the error box was then performed 
locally varying the mass, initial rotational velocity, and age of the models. The mass was 
varied around 7.10 M0, in particular from 7 to 7.20 M©, in steps of 0.01 M©. Rotation 
was considered under two main assumptions: uniform rotation and differential rotation (see 
Section [2] for more details). For both assumptions, the rotational velocities considered range 
from to 5 to 20kms~^ at the stellar surface. Then, the modelling was refined by analysing 
the rotational splittings and their asymmetries. Rotational splittings are defined as 

S = ^(^^+1 - z^-i) (3) 

and asymmetries of these splittings are defined as 

A = iy_i + - 2z/o . (4) 

where sub-indices ±1 represent the value of the azimuthal order m for a given i. The 
information provided by the asymmetries can be completed by the semi-splittings, which 
corresponds with A+ = z/+i — uq, and A~ = z/q — z/_i, respectively. Both quantities are 
related such as Eq. H] transforms into 

A = A+-A-. (5) 

Moreover, this analysis also permits to extract information about the internal rotation profile 
of the star, since, as we discuss in the next sections, S and A are both sensitive to changes 
in the rotation profile. 



-17- 



5. Results and discussion 

Models with a mean density of about p = 0.064 gcm~^ were found to fit better the 
observed frequencies. The models within ±0.05 M0 of the mass value calibrated using non- 
rotating models yield similar results, but for slightly different ages. Due to the so small 
rotational velocity, it is plausible to assume the same omlt and rfov parameters which were 
calibrated by the non-rotating models. Therefore, only the mass, age, and initial rotational 
velocity of the models were varied. 

Models were then selected to fit at least the observed frequencies /i (identified as the 
fundamental radial mode) and the triplet (/s, /4, /2), identified as gi. With these criteria, a 
model with a mass of 7.13 Mq, a rotational velocity (in the surface) of about 7kms~^, and 
an age of about 14.9 Myr (which corresponds with logTefj = 4.365) was found to better fit 
the observational frequencies. As expected, when rotation is taken into account, the stellar 
parameters of the models are similar to those of the non-rotating best model (NRl). For 
each type of rotation, i.e. uniform and shellular rotation, the models better matching the 
observed frequencies, splittings, and asymmetries, are URl and SRI, whose characteristics 
are summarised in Table [2] (their corresponding list of frequencies are reported in Table [3]) , 
and whose internal rotation profiles are depicted in Fig. |4l Figure [5] shows the 
evolution of the theoretical frequencies of the gi, pi and p2 ^ = 1 triplets for the selected 
model (assuming shellular rotation). 

With these models, nine of the observed frequencies were identified as two i = 1 triplets: 
a gi triplet composed by the observed frequencies (/a, f^, /2), a pi triplet which corresponds 
with {fi2, fe, fj), and a p2 triplet which corresponds with (/5,/9,/io)- Furthermore, fi was 
identified as the fundamental radial mode. For the remaining frequencies (/§ and 
/ii), the selected models present similar predictions. While /§ is identified as a 
(n = 0,£ = 2), all the models match /n with a i >7 mode. In particular the best 
identification found corresponds to {n,£,m) = (—3,9,8). 

In general, these results are compatible with previous studies that can be 
found in the literature, except that we do identify and use (for asteroseismic 
purposes) all the non -axisymmetric components of t he £ = 1 triplets, in particular 



P2. Differences with IDziembowski &: Pamyatnykhl (120081 ) 's work are principally 



related to this last split mode, whose identification depends strongly on the use 
of standard or non-standard models. 
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Table 3. List of theoretical frequencies (d~^) for the three selected best models NR, SRI, 
SR2, and UR (see Table [2]). The model A does not take rotation into account, hence only 
the frequencies of the m = components are reported. 



Mode z/NR,i z/sRi,i z/sR2,i jyvR,i 



FO 

gi,-i 

gl,0 

gi,+i 
Pi,-i 
Pi,o 
Pi,+i 

P2-1 

P2,0 

P2,+l 



5.75818 
5.63699 

6.24442 

7.88829 



5.74396 
5.65407 
5.63402 
5.61397 
6.26192 
6.24026 
6.21785 
7.89759 
7.87398 
7.84899 



5.72332 
5.70359 
5.67979 
5.65585 
6.31067 
6.28685 
6.26212 
7.87980 
7.85263 
7.82361 



5.77681 
5.66402 
5.64512 
5.62652 
6.27800 
6.25494 
6.23102 
7.94832 
7.91898 
7.88731 
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5.1. Analysis of rotational splittings and asymmetries 

In Fig. [6l the predicted semi-splittings (left column) and asymmetries (right column), 
are compared with the corresponding observed values given in Table HI For all panels in that 
figure, the shaded vertical regions indicate the range of Teg values given by the best models 
SRI and URl (see Table [2]). 

This comparison is performed simultaneously for models evolved assuming uniform and 
shellular rotation profiles (hereafter UR and SR models). For the sake of brevity, in the fol- 
lowing, A^, A°, are used to represent the predicted and observed asymmetries, x representing 
the triplets gi, pi, or p2. Similarly, A+'^ and A+'~ are used to represent the predicted and 
observed semi-splittings. Let us now examine each selected triplet separately. 



5.1.1. The gi triplet 

Analysis of the gi triplet (Fig. [6l left top panel) reveals that for the SR models, the 
predicted lA+l remain lower than |A~| for almost the whole range of effective temperature 
studied (almost identical in the shaded region). On the other hand, for UR models, |A~| are 
predicted lower than |A+| for effective temperatures higher than 23014 K (log Teg = 4.362), 
approximately. The deviation of the semi-splitting predictions (in the shaded region), defined 
as eA+,- = A+~ - A+~, is about 2 x 10"^ d~^ for the SR models and 3-4 x 10"^ d~^ for 
the UR models. On the other hand, cooler SR models would fit better the observations 
(eA+.- = 10~^d~^ for models around log Teg = 4.35). In general, the evolution of the semi- 
splittings for both types of rotation is found to be quite different. This is somehow expected 
since g modes are very sensi t ive to variations of the rotational velocity of the core. In fact, it 



was shown by ISuarez et al.l (120061 ) that a shellular rotation profile modifies significantly the 
radial displacement eigenfunctions, especially for g and mixed modes. This implies that the 
presence of shellular rotation may affect both the rotational splitting itself and its asymmetry. 
The general definition of the first-order rotational splitting kernel can be written as 

(2 2/01^0 + zl) + r/o(2/oi + A^o - 2 VmZQ - zl) 



4 



R 

2 



(6) 







where A = and y and z represent the vertical and horizontal displacement normalised 

eigenfunctions, respectively. The second term within the square brackets accounts for the 
presence of shellular rotation through r^o, defined in Eq. [H When a uniform rotation is 
considered, this term becomes null. In Fig. [7| such kernels are depicted for the gi, pi, and 
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P2 triplets. Notice that, when considering a shellular rotation profile, a bump in the energy 
distribution near the yU-gradient zone (see Fig. Hj) comes up, for gi, at the expense of the 
energy of the outer layers, which could explain the different behaviour of the semi-splittings 
predicted by the UR and SR models. 
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Table 4. Observed semi-splittings, A^), and asymmetries, A" of the i = 1 triplets gi, 
Pi and p2, as identified using a model of 7.13 which rotates with v ~ 5.5kms~^ in the 
surface. Quantities are giv en in d~^. Uncertainties, S , are calculated from data in 



Jerzykiewicz. et al.l (120051 ). 



Mode 


A+ 


A- 






<5(A+) 


6{A~) 


6{A^) 


gl 


0.016630 


-0.01723 


-0.000598 




±6 X 10-6 


±7 X 10-^ 


±1 X 10"^ 


Pi 


0.0191 


-0.0203 


-0.0012 




±2 X 10-^ 


±2 X 10~^ 


±2 X 10"^ 


P2 


0.0161 


-0.0156 


+0.0005 




±1 X 10"^ 


±2 X 10"^ 


±2 X 10"^ 
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Concerning the asymmetries, A^^ are predicted to diminish for decreasing effective tem- 
perature (Fig. [6l right panel on top). In other words, asymmetries decrease while the star 
evolves. For UR models, such a decreasing is more rapid than for the SR models, and the 
asymmetries fit the observed value at logTefj ~ 4.357, which represents a difference of ~ 400 
K with respect to the models that fit the observed frequencies (shaded area). On the other 
hand, the asymmetries predicted by the SR models never fit exactly the observed value. In 
the shaded region, the deviation of the asymmetry defined as = — is found to 
be 6 X 10~^d~^ for the SR models whereas for UR models eA = 10^^ d^^. Such deviations 
are respectively, two and three orders of magnitude larger than the observed asymmetry 
uncertainty. 

As for the gi semi-splittings, better results are found for cooler models, especially for 
UR models around logTgfj = 4.357, and for SR models near log Teg = 4.35. In any case, 
for low order g and p modes, asym metries are sensitive to variations of the rotation profile 



near the core. Indeed, as shown by ISuarez et al.l (120061 ). the analytical form of A (using a 



perturbative theory) can be written as the second-order term 

4A — 6uJq 

where VL is the rotational velocity at the stellar surface, and uoq is the unperturbed oscillation 
frequency. The Jc integral contains a complex combination of structure and oscillation 
terms which are modified by the rotation profile and its derivatives. Analysis of this term 
is definitely necessary to construct simplified kernels for A (work in progress). Such kernels 
will help us to better understand the behaviour of the asymmetries, and, especially their 
sensitivity to variations of the internal rotation profile. 
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Fig. 4. — Rotation profile normafised to the rotation velocity in the surface, fig, as a function 
of the normalised stellar radius, for the best SR model. For comparison, the rotation profile 
of UR model with similar is also depicted. For illustration purposes, the Brunt- Vaisala 
frequency for the SR model (re-scaled as to fit the plot) has been used to indicate the location 
of the //-gradient zone (r ~ 0.15 — 0.2 R0). 



5.1.2. The pi triplet 



For pi, the predicted curves of |A+| remain lower than those of |A^| in the whole 
range of effective temperature studied, which is compatible with the observations (Fig. [6l 
left middle panel). Notice that, in the shaded region, the SR models predict absolute values 
for the semi-splittings closer to the observed ones {e^+.- ~ 10""^ d~^) than those predicted 
by the UR models, for which eA+.- = 4 x 10^^ d^^. However, this situation is reversed for 
the splitting asymmetries (Fig. [6l right middle panel). In particular, around logTes = 4.365, 
the asymmetries predicted by the UR models {eA = 6 x 10~^d~^) are slightly closer to the 
observed values than those predicted by the SR models (e^ = 2 x 10~^d~^). Even so, such 
differences are of the order of magnitude of the observational uncertainty of A^-^ , which makes 
it difficult to discriminate between both types of rotation. Contrary to the gi results, cooler 
models do not fit better the observations. 



5.1.3. The p2 triplet 

In the case of p2, the predicted semi-splittings are larger than the observed value in the 
whole range of effective temperatures, except for log Teg ~ 4.35, for which SR predictions 
are almost coincident with the observations (Fig. [6l left bottom panel). In the shaded 
region eA+,- = 6 x 10"^ d~^ for the SR models and eA+,- = 1-6 x 10~^d~^ for the UR 
models. Similarly to the gi case, the best results are given by the SR models for effective 
temperatures around the cooler limit of the photometric uncertainty box (log Teg = 4.35). 



4.365 4.360 4.355 4.350 



Fig. 5. — Evolution of the theoretical oscillation frequencies corresponding, from to bottom, 
to the £ = 1 triplets p2, pi and gi, for the selected 7.13M0 model. Horizontal lines represent 
the observed frequencies, bottom to top, f^, /g, and /g, identified as the m = components 
of the triplets p2, Pi and gi, respectively. 
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Furthermore, it is worth noting that contrary to the pi and gi cases, the observed semi- 
sphtting |A~| is smaller than |A+|, whereas both SR and UR models predict |A^| > |A+| 
(Fig. [HI left bottom panel). This results in a positive observed asymmetry, whereas both 
UR and SR models predict negative values (Fig. [HI right bottom panel). In the vicinity of 
logTeff = 4.365 (shaded region), the difference between the observed and predicted asymme- 
tries is about eA = 10"^ d"^ (e^ = 1.5 x 10"^ d"^ for the SR models and = 2.8 x 10"^ d"^ 
for the UR models), which represents a difference of one order of magnitude with respect 
to the observed value. Such apparently marginal contradiction between predictions and ob- 
servations could be a consequence either of 1) an incorrect mode identification, that is, the 
observed frequencies concerned do not belong to the rotationally split mode, or 2) the use 
of a wrong description for the rotation profile, particularly in the outer shells of the star 
(note the small influence of the selected rotation profile near the /^-gradient zone. Fig. [7]). 
A priori, none of the these possibilities can be discarded. In order to solve this problem, 
improvements on both the observations and modelling are required. From the observational 
side, an improvement of accuracy with which the concerned observed frequencies are deter- 
mined might help to confirm the observed asymmetry and its sign. Moreover, the detection 
of additional frequencies (e.g. with the help of space missions), may provide new insight on 
the current mode identification. 

From the theoretical side, the second possibility given above is related with angular 
momentum redistribution, which plays an important role. In particular, balance between 
rotationally induced turbulence and mer idional circ ulation generates mixing of chemicals 



and redistribution of angular momentum (jZahrull992l ). which affects the rotation profile and 



the evolution of the star. The present technique is, therefore, especially suitable for testing 
that theory by providing estimates for the coefficients of turbulence using only asteroseismic 
observables. This can be illustrated by artificially modifying the physical conditions beyond 
the convective core, which can be done varying the overshooting parameter. 
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Fig. 6. — Evolution of the theoretical semi-splittings (left column) and asymmetries (right 
column) as a function of the effective temperature (logarithmic scale) for the selected 7.13 Mq 
model with a rotational velocity of 7kms"^ , approximately. The observed values are de- 
picted with horizontal lines, and the corresponding observational uncertainties are repre- 
sented by horizontal shaded bands. Each panel row shows the results, from top to bottom, 
for the £ = 1 triplets gi, pi, and p2, respectively. In left panels, solid hues represent the 
positive semi-splittings A"*", and dashed lines, the negative ones A~. For the asymmetries, 
the results obtained from UR and SR models are represented by dashed and solid lines, 
respectively. The shaded vertical region indicates the range of effective temperatures defined 
by the models SRI and URl (more details in Section [5?II) . 
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As expected, such variations principally affect the low-order gi and pi (see Fig. [7] for 
a comparison between the rotational kernels of UR, SR-ciov = 0.28, and SR-ciov = 0.24 
models). In Fig.[H], the results for the asymmetries given in Fig. [H] (right column), which were 
computed with an overshooting parameter of d^^ = 0.28, are compared with those obtained 
from SR models computed with dov = 0.24. The best model found for this overshooting 
parameter value is SR2 (see Tables [2] and [3]) . In particular, a variation in dov of 0.04 results 
in differences in the asymmetry of the order of those found between UR and SR {d^v = 0.28) 
models, i.e. a few 10~^d~^ for gi and pi, and about 10~^d~^, which represents almost one 
order of magnitude smaller than the difference between the results yield by the UR and SR 
((iov = 0.28) models. Moreover, in the case of pi, the SR {d^v = 0.24) models fit the observed 
asymmetry in the limits of the observational uncertainties. 
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r/R 



Fig. 7. — First-order rotational splitting kernel, from top to bottom, for the gi, pi, and 
P2 modes for the selected uniformly-rotating models (UR, dashed line) and differentially- 
rotating models (SR). For the latter models, two dov values are considered: 0.28 (continuous 
hne) and 0.24 (dot-dashed line). Note that the position of the zeros and maxima for UR 
models is shghtly shifted with respect to those of the SR models. This can be explained by 
small differences in the radii of both models, which, as expected, principally affects to the 
gi and pi split modes. 
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Fig. 8. — Comparison of the splitting asymmetries shown in Fig. [5] (for models with dow = 
0.28), with those given by the models computed with dow = 0.24 (dot-dashed line). The 
additional vertical shaded band represents the range of T^a defined by the models SRI, SR2, 
and URl. 
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Furthermore, recent theoretical studies (Andrade, Suarez & Goupil 2008, work in progress), 
seem to indicate that even rough variations of the rotation profile may modify the asymme- 
tries of the split modes. This modification can be as important as to change the sign of the 
asymmetries. According to this, the frequencies around p2 cannot be discarded as belonging 
to a rotationally split mode. However, if confirmed, the presence of such different asymme- 
tries in the oscillation spectrum is very interesting, since they may significantly constrain 
the models providing information about the structure and rotation profile in the zones where 
they have more amplitude. 

6. Conclusions 

An astcroscismic analysis of the f3 Ccphci star u Eridani is presented, with focus on 
the study of the internal rotation profile. To this aim, a new method is presented based 
on the analysis of rotational splittings and their asymmetries. Some of the most updated 
asteroseismic modelling techniques are used, in particular, analysis of mode stability and 
improved descriptions for rotation effects. Regarding the latter, the so-called pseudo-rotating 
models are used, which consider radial differential rotation profiles (shelluUar rotation) in 
both the evolutionary models and adiabatic oscillations computations. This represents an 
important qualitative step with respect to previous theoretical works. 

The present work is divided into two parts. In the first part, a comparison of the differ- 
ent numerical packages (Liege, Warsaw-New Jersey, and Granada packages) was performed 
(the latter has been used in the present work). This comparison was performed using a 
mass-metallicity relation of the models fitting two of the observed frequencies. In that case, 
important differences in mass, up to 0.25 Mq between Liege models and ours, were found. 
Interestingly, differences between ASTA04's models and ours can be explained by the dif- 
ferent treatments of the overshooting implemented in the evolutionary codes. ASTA04's 
models only take into account density variations in the overshooted region, which slightly 
affects the temperature gradient. On the other hand, our models were constructed consid- 
ering an additional restriction by imposing that the real temperature gradient must equal 
the adiabatic gradient. Physically, this implies that, either the heat is transported efficiently 
outwards from the stellar core through convective movements due to the overshooting (the 
latter case), or purely by radiation (ASTA04's treatment). This constitutes a very interest- 
ing challenge for astcroscismology, since the oscillation modes (low-order g and p modes) are 
sensitive to the physical description of the /i gradient zone. 

Then, the previous exercise is extended to four frequencies with mode excitation in- 
cluded. In that case, models were also constrained to match the fundamental radial mode 
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and three m = splitting components. Non-standard models built with a significant de- 
crease of the initial hydrogen abundance, = 0.50, were necessary to match and excite the 
modes. This provided models similar, but not to those found by ASTA04 using 

a solar iron relative abundance, and the same initial hydrogen abundance. The remaining 
differences between ASTA04's models and ours indicate that, for this level of the model 
accuracy, the present modelling still depend on the core overshooting. This parameter may 
change the physical conditions beyond the convective core. In that region, other physical 
processes take place, such as the mixing of chemical elements due to rotation. In particu- 
lar, balance between rotationally induced turbulence and me ridional circ ulation generates 



mixing of chemicals and redistribution of angular momentum (jZahrull992l ). This affects the 
rotation profile and the evolution of the star. 

Secondly, the method here presented studies the asymmetries of the split modes in order 
to refine the modelling, and more importantly, provide information about the rotation pro- 
file of the star. To do so, pseudo-rotating models were built using the physical parameters 
provided in the first part of the work. Models with masses around 7.13 Mq, and ages around 
14.9 Myr, were found to fit better 10 of the 14 observed frequencies, which were identified 
as the fundamental radial mode and the three i = 1 triplets gi, pi, and p2. For these 
modes, a comparison between the observed and predicted splittings and their asymmetries 
was performed. Two type of rotation profiles were considered: uniform rotation and shel- 
lular rotation profiles. Differences between predictions and observations were found to be 
of the order of 5 x 10"'^ and 10~^d~^, for the rotational splittings and their asymmetries, 
respectively. For this last mode, none of the selected models reproduce neither the splittings, 
generally larger than the observed ones, nor the asymmetries, whose predictions have the 
opposite sign than the observed values. Although this result might indicate that the frequen- 
cies around p2,'m = mode do not belong to this rotationally-split mode, other possibilities 
cannot be discarded. In fact, a wrong physical description of the rotation profile may be 
responsible for a so peculiar result, particularly in the outer regions of the star. This result is 
very important because, up to now, none of the physical phenomena in which rotation plays 
a role predict variations of the rotation profile in that region. Furthermore, it is shown that 
asymmetries are quite dependent on the overshooting of the convective core. Therefore, the 
method here presented is suitable for testing the theories describing the angular momentum 
redistribution and chemical mixing due to rotationally-induced turbulence. 

In general, the seismic models which include a description for shellular rotation yield 
slightly better results as compared with those given by uniformly-rotating models. Even so. 



^Cf. the explanation given above when matching only two frequencies 
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further improvements are necessary in order to better constrain the modelling of v Eridani. 

In particular, efforts should be focused in searching a better description of the rotation 
profile, for which a detailed study of asymmetries is required (work in progress). This may 
be enhanced by analysing other (5 Cephei stars with larger rotational velocities, so that 
splittings and asymmetries are not of the same order than the frequency uncertainties. 
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